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Abstract: Vegetation roots contribute to soil fixation and reinforcement, thus improving soil resistance 
against erosion. Generally, the amount of soil fixation presented by roots mainly depends on root density 
and tensile strength. In the present study, we conducted the research in order to further understand the 
biotechnical properties of Haloxylon persicum and also to quantify its role in increasing soil cohesion in arid 
lands of Iran. Ten H. persicum shrubs were randomly selected for root distribution and strength 
investigations, in which five samples were set on flat terrain and other five samples on a moderate slope 
terrain. The profile trench method was used to assess the root area ratio (RAR) as the index of root 
density and distribution. Two profiles were dug around each sample, up and downslope for sloped 
treatment and north and south sides for flat treatment. The results showed that RAR increased with 
increasing soil depth and significantly decreased in 40—50 cm layers of downhill (0.320%) and 50—60 cm 
for uphill (0.210%). The minimum values for the northward and southward profiles were 0.003% and 
0.003%, respectively, while the maximum values were 0.260% and 0.040%, respectively. The relationship 
between the diameter of root samples and root tensile strength followed a negative power function, but 
tensile force increased with increasing root diameter following a positive power function. The pattern of 
increased cohesion changes in soil profile was relatively similar to RAR curves. The maximum increased 
cohesion due to the presence of roots in uphill and downhill sides were 0.470 and 1.400 kPa, respectively. 
In the flat treatment, the maximum increased cohesions were 0.570 and 0.610 kPa in northward and 
southward profiles, respectively. The analysis of variance showed that wind and slope induced stresses did 
not have any significant effect on the amount of increased cohesion of H. persicum. The findings served to 
develop knowledge about biotechnical properties of H. persicum root system that can assist in assessing the 
efficiency of afforestation and restoration measures for erosion control in arid lands. 
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1 Introduction 


Vegetation plays an important role in slope stabilization and erosion control through mechanical 
and hydrological effects (Genet et al., 2008; Mao et al., 2012). Among different parts of 
vegetation, roots are reported to have an essential influence on slope stabilization and erosion 
control (Wu et al., 1979), mainly due to their mechanical effects (Simon et al., 2006). Roots can 
improve some engineering characteristics of soil, such as cohesion and shear strength, which is 
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known as soil fixation (Genet et al., 2006) or root reinforcement (Genet et al., 2008). The 
presence of roots as a material with strong tension in soil with weak tension provides a composite 
material with enhanced strength (Pollen, 2007; De Baets et al., 2008) and higher shear strength. 
The higher soil shear strength will result in a higher resistance against erosion and instability. 
Many efforts have been made to quantify the effect of roots in changes of soil cohesion and some 
models have been developed for this purpose, such as the model of Wu et al. (1979), the RipRoot 
model (Pollen and Simon, 2005) and the Root Bundle model (Schwarz et al., 2013). Among them, 
the model of Wu et al. (1979) is the simplest and requires minimal input parameters (Schwarz et 
al., 2012). This model estimates the increased cohesion as a result of vegetation's root presented 
in the soil using root density and strength parameters (Greenway, 1987; Gray and Sotir, 1996). 
Previous researches concerning soil bioengineering reported root density and strength as the main 
biotechnical characteristics of different vegetation covers that control the fixation or 
reinforcement effect of roots (Gray and Sotir, 1996). 

Root area ratio (RAR) has been used as the most popular indicator of root distribution and 
density by many researchers (e.g., Abernethy and Rutherfurd, 2000; Genet et al., 2008). Also 
tensile strength is considered as an important factor that governs the amount of soil fixation and 
stabilization and, therefore, many researchers have investigated this parameter (Bischetti et al., 
2005; Pollen, 2007; Genet et al., 2008). Both of the above mentioned parameters are affected by 
many biotic and abiotic factors, such as plant species, ages, and densities (Gray and Sotir, 1996; 
Nilaweera and Nutalaya, 1999; Schmid and Kazda, 2001; Genet et al., 2008). Knowledge of root 
properties, such as density and strength, is important for the implementation of soil 
bioengineering measures and corresponding efficiency assessment, and can be used in quantifying 
the effect of different plants in erosion control measures (Pollen, 2007). 

In recent years, many studies have been conducted regarding soil bioengineering in different 
ecosystems, but there is still a lack of information about some ecosystems, especially arid 
ecosystems and corresponding plant species. Some studies in recent years have showed a higher 
level of resistance to erosion of soils that are covered with vegetation (Gray and Sotir, 1996; De 
Baets et al., 2008), therefore, soil bioengineering techniques can be used for erosion control (Gray 
and Sotir, 1996), and soil fixation (Genet et al., 2006) in arid lands. 

The area of deserts in Iran is about 21x10° hm? and nearly 70% of the country is covered with 
arid and semi-arid lands (Khosroshahi et al., 2011). Over the past century, extensive areas of arid 
and semi-arid regions (including Iran) have been gradually converted to afforestation (Jazirei, 
2009). H. persicum (white saxaul) is one of the most predominantly selected broadleaved shrubs 
for afforestation/reforestation in the arid and semi-arid regions of Iran, Iraq, Afghanistan, and 
China (Pyankov et al., 1999; Jazirei, 2009; Nosrati et al., 2013). This is because it tolerates 
droughts and high temperature extremes (Jazirei, 2009). Also, this species is conventionally used 
as livestock feed and firewood, and is important for stabilizing sand dunes (Tobe et al., 2000; Dai 
et al., 2015). In the current study, we aim to further understand the biotechnical properties of H. 
persicum and also quantify its role in increasing soil cohesion that may lead to soil fixation and 
increasing soil resistance against erosion. 


2 Materials and methods 


2.1 Study area 


The study region with an area of 85x10° hm? is located in a small town near Zavareh in central 
Iran. Average elevation is 1125 m a.s.l. and the slope ranges from 0%—35%. Prevailing wind 
direction is south and the amount of annual precipitation is 110 mm. The regional climate is hot 
and dry and the temperature differences between the day and night are high. Shrub spacing of H. 
persicum plantations is about 3 mx3 m (1100 stems/hm?). The average collar diameter of samples 
was 14 cm. 


2.2 Root distribution assessment 


Ten H. persicum shrubs were randomly selected in a 43-year old plantation of H. persicum, for 
root biotechnical properties investigations. The profile trench method (Böhm, 1979; Bischetti et 
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al., 2005; Abdi et al., 2010) was used to assess the root distribution of H. persicum shrubs. In this 
study, RAR, the ratio of the sum of the root areas to the area of soil profile that they intersect, was 
chosen as the index of root distribution in the soil profiles. To assess the effect of slope on the 
biotechnical properties of the root system, we set five shrub samples on a flat terrain and other 
five samples on a south-facing slope with a gradient of 35%. Two profiles were dug around each 
sample shrub. In the sloped treatment, one profile was dug upslope side (northward) and the other 
was dug downslope side (southward) of each shrub sample. Also on the flat treatment, profiles 
were dug on the north and south sides of each sample (Fig. 1). As the prevailing wind direction in 
the study area is southward and considering shrub samples as reference points, northward and 
southward profiles were considered as windward and leeward, respectively. Therefore, the effect 
of wind and slope on the biotechnical properties of the samples can be assessed. 


(a) Flat terrain (b) Slope terrain 
Northward e nee 
profile Wind direction Uphill 
E ELF Loriotile 
ae leshrab #8Sample a. sl 
ee Northward shrub 
Southward Dow nhi 
profile profile 


Fig.1 Schematic design of sampling on flat (a) and slope terrains (b) 


The profiles were excavated at a distance of 20 cm from the collar of samples and profile 
excavations were continued until no roots were seen on trench walls. The horizontal dimensions 
of the profiles were approximately 100 cmx50 cm. Layers of 10 cm depth were marked on the 
vertical trench walls near to samples and all root diameters were measured by a digital caliper. 
Using the following equation, the area occupied by roots in each soil layer (i.e., 4, in model of Wu 
et al. (1979)) was determined: 


AZ Exa), (1) 


where 4, is the sum of root areas (mm/7); and d; is the diameter of the i" root (mm). 
RAR values were calculated by dividing the sum of the root areas in each 10 cm horizon (4,) to 
the soil area (A), i.e., 1000 mmx100 mm in this study. 


2.3 Root strength tests 


Live roots for mechanical tests were randomly collected from trenches on both sides of the 
sample shrubs at a depth of about 30 cm as suggested by Cofie and Koolen (2001). To preserve 
samples from fungal contamination, we used the treatment suggested by Bischetti et al. (2005). 
Therefore, the root samples were sprayed with a 15% alcohol solution and samples were stored in 
plastic bags. The tensile tests were carried out on fresh roots within three days after the sampling. 
Suitable lengths of root samples were cut to about 150 mm lengths for strength tests. The 
diameter of the root samples was measured in three different positions along the root length and 
the average was considered as the root diameter. The mechanical properties tests were carried out 
using an Instron Universal Testing Machine (Model 4486, UK). The test speed was 10 mm/min 
(Bischetti et al., 2005; Mattia et al., 2005; Pollen, 2007). The root ends were clamped to jaws and 
the tests were continued until rupture occurred in the root samples. Tensile strength was 
calculated by dividing the maximum force required to rupture the root samples (N) by the root 
cross-sectional area (mm°): 

T, ===, (2) 


Tyd? 
4 


where T, is the tensile strength (MPa); Fmax is the maximum force required to rupture the root 


Ehsan ABDI et al.: Soil fixation and erosion control by Haloxylon persicum roots in arid lands, Iran 


samples; and d is the root diameter (mm). 


2.4 Soil strength parameters 

To estimate the shear strength parameters of soil for the study area, we conducted a direct shear 
test on 50 mmx50 mm soil samples. The direct shear test was performed with three replications 
with normal stresses of 10, 20 and 30 kPa. The test speed was fixed to 0.5 mm/min. Based on the 
results, the values of C and ọ were 0.08 kPa and 22°, respectively. 


2.5 Increased soil cohesion 


The model of Wu et al. (1979) was used to estimate the potential increased soil cohesion (C,) due 
to the presence of roots. This increased cohesion can result in an increase in soil shear strength. 
The base of the root reinforcement model of Wu et al. (1979) is the Coulomb equation that 
calculates soil shear strength based on two main parameters, cohesive and internal friction angle 
of soil, soil shear strength can be calculated as follows: 

S=C+0, tang, (3) 
where S is the soil shear strength (kPa); C is the soil cohesion (kPa); on is the normal stress on the 
shear plane (kPa); and ø is the soil internal friction angle (°). 

In the model of Wu et al. (1979), the modified Coulomb equation for root-permeated soils 
becomes: 
S=C+C, +0, tang, (4) 
where C, is the increased cohesion due to the presence of roots (kPa). 
C, can be represented by: 
C. =kxt,, (5) 
where Ń is a constant factor; and ¢, is the mobilized root tensile strength per unit area of soil. 
The & factor can be represented by the following equation: 
k =sin@+ cos Otan ø, (6) 


where @ is the angle of shear distortion in the shear zone (°); and ø is the soil internal friction 
angle (°). 

In the current study, the angle of shear distortion is assumed to be equal to 45° regarding the 
result of Waldron (1977) and the friction angle of the site soil was 22°, then, k factor in Equation 
6 is taken as 1.03. 

The mobilized root tensile strength per soil unit area (¢-) can be written as: 

t_=T.Xa,, (7) 
where T, is the average tensile strength per average root cross-sectional area (MPa); and a; is the 
RAR computed as A,/A, where A, is the total cross-sectional area of all roots (mm7*), and A is the 
area of soil (mm7’). 

As tensile strength and root diameter have a reverse power relation and considering the 
variability in the diameter of roots, we changed Equation 7 as follows: 


N 
c= 37, (8) 
m A 
where i is the diameter class; and N is the number of classes. 
In the current study, by combining the root strength measurements and the RAR values 
calculated from trenches beside H. persicum samples, we calculated the potential increased 
cohesion of soil due to the presence of H. persicum roots using Equation 8. 


2.6 Statistical analysis 


SPSS 16.0 statistical software was used to analyze the data. The Kolmogorov—Smimov test was 
used to check the normality of the data before proceeding with the analyses and as the normality 
assumption was violated (P<0.05) and data were transformed. Paired samples T test was used to 
compare RAR and C, values of two side profiles. Analysis of covariance (ANCOVA) and analysis 
of variance (ANOVA) were used to investigate the differences in live root tensile strength (with 
root diameter sets as the covariate factor) and increased cohesion served by roots among different 
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treatments, respectively. Duncan multiple range test was used to compare mean increased 
cohesion among different treatments. 


3 Results 


3.1 Root distribution 


As standard error bars shown in Figure 2, the RAR values showed a great variability considering 
soil depth and hillside. Generally, RAR values increased with increasing soil depth and then 
decreased again, with maximum values being in the middle soil layers (Fig. 2). The minimum 
values of RAR were in the 0-10 (0.003%) and 80—90 cm depth (0.006%) for uphill and downhill, 
respectively. The maximum values were in the 40-50 cm layer for downhill (0.320%) and 50-60 
cm layer for uphill (0.210%). The differences of RAR values were only significant in the fourth 
layer of soil with higher values belonging to downhill profiles. 


RAR 
0.0 0.1 0.2 0.3 0.4 0.5 0.6 


1l0yaa --e-- Downhill profile 
204 aci —e— Uphill profile 


Depth (cm) 
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Fig. 2 Relationship between root area ratio (RAR) and soil depth for downhill and uphill profiles. Bars indicate 
standard errors. Means with the same lowercase letters are not significantly different between soil profiles 
(P>0.05). 


The values of RAR in flat treatment in both northward and southward profiles showed smaller 
values and less variability compared with sloped profiles, except the last layer of northward 
profiles. The RAR values was increasing to middle depth and then decreasing to deeper layers 
with an exception in the last layer of northward profiles. The minimum values for the northward 
and southward profiles were 0.003% and 0.003%, respectively, while its maximum values were 
0.260% and 0.040%, respectively (Fig. 3). The maximum values for northward and southward 
occurred in the eleventh layer (90-100 cm) and the third layer (20-30 cm), respectively. The 
difference of RAR values was only significant in the eleventh layer of soil, with higher values for 
northward profiles. 


0.0 0.1 0.2 0.3 0.4 0.5 
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a —e— Northward profile 
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Fig. 3 Relationship between root area ratio (RAR) and soil depth for southward and northward profiles. Bars 
indicate standard errors. Means with the same lowercase letters are not significantly different between profiles 
(P>0.05). 
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3.2 Root strength tests 


In total, 91 root specimens were tested for different treatments. The root diameter for sloped 
treatment varied between 1.75-8.61 mm for uphill and 1.59-6.57 mm for downhill, the mean 
tensile force values were 17.02 (+1.83) N (uphill) and 17.52 (+1.56) N (downhill), the minimum 
and maximum recorded values were 2.29 and 40.07 N for uphill and 3.67 and 41.58 N for 
downhill, respectively (Fig. 4). The tensile force increased with increasing root diameter 
following a power function as shown in Figures 4 and 5. 
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Fig.4 Relationship between tensile force at rupture point and root diameter for downhill and uphill profiles 


The range of the root diameter of flat treatment varied between 1.52 and 5.62 mm for 
northward profile and 1.61-5.48 mm for southward profile, the mean tensile force values were 
12.93 (+1.79) N (northward profile) and 12.05 (+2.49) N (southward profile), and the minimum 
and maximum recorded values were respectively 4.15 and 32.93 N for northward profile and 3.71 
and 43.61 N for southward profile (Fig. 5). 
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Fig. 5 Relationship between tensile force at rupture point and root diameter for southward and northward 
profiles 


The power function was fitted to scatter plot of tensile strengths (y axis) versus root diameter (x 
axis). The results showed that power function was statistically significant in all treatments 
(P<0.05). The relation followed a negative power function and tensile strength increased with 
decreasing root diameter (Figs. 6 and 7). The ranges of tensile strength were 0.15-6.31 MPa for 
uphill, 0.56-5.56 MPa for downhill in sloped treatment, 0.61-4.92 MPa for northward and 
0.62-5.49 MPa for southward in flat treatment. 

The coefficients of power function parameters for different treatments and statistical 
information are shown in Table 1. The results of ANCOVA showed that there is no significant 
differences between the tensile strength of different treatments (P>0.05) and only the effect of the 
root diameter as a covariate factor was significant (P<0.01). 
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Fig.6 Relationship between root tensile strength and root diameter for uphill and downhill profiles 
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Fig. 7 Relationship between root tensile strength and root diameter for southward and northward profiles 


Table 1 Model summary of power functions fitted to tensile strength and root diameter 


Profile Number of samples a B R? P value 

Uphill 31 6.02 -1.30 0.416 0.000 
Downhill 29 5.25 0.97 0.404 0.000 
Northward 15 7.55 -1.36 0.492 0.000 
Southward 16 3.36 -0.75 0.130 0.000 


3.3 Increased soil cohesion 


The pattern of increased soil cohesion (C,) was relatively similar to the pattern of RAR. 
Regarding sloped treatment, the ranges of increased cohesion in uphill and downhill sides were 
0.020—0.470 and 0.005-1.400 kPa, respectively (Fig. 8). The maximum increased soil cohesion 
was presented in the fifth layer and the differences of C, values were significant in the third, 
fourth and seventh layers of soil, with higher values for downhill profiles. 

In flat treatment, in northward profiles, cohesion rate increased to 0.570 in the sixth layer. The 
soil reinforcement in the southward profile was the maximum in the third layer with 0.610 kPa 
(Fig. 9). The differences of C, values were only significant in the seventh layer of soil with higher 
values for northward profiles. 

The results of ANOVA showed that there were significant differences in increased cohesion 
among treatments (F=3.87, P<0.05). The results of Duncan multiple range test showed that 
downhill profiles had significantly higher values than the other profiles (Fig. 10), but the 
differences among uphill, leeward and windward were not statistically significant (P>0.05). 
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Fig. 8 Increased cohesion (C,) due to presence of roots for uphill and downhill profiles under different soil 
depths. Bars indicate standard errors. Means with the same lowercase letters are not significantly different 
between soil profiles (P>0.05). 


Depth (cm) 


Cer aAUsWHY 
ooocoocoooce 


--@-- Southward profile 
—e— Northward profile 


Fig. 9 Increased cohesion (C+) due to presence of roots for northward and southward profiles under different soil 
depths. Bars indicate standard errors. Means with the same lowercase letters are not significantly different 
between soil profiles (P>0.05). 
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Fig. 10 Mean increased cohesion (C+) under different treatments. Means with the same lowercase letters are not 
significantly different among different treatments (P>0.05). 


4 Discussion 


Many biotic and abiotic factors can affect root properties. Forest management, genetic factor, soil 
and climate characteristics (Bischetti et al., 2005), age and stand density (Genet et al., 2006) are 
reported to affect root density and distribution. In the current study, high variability of RAR 
values was observed. Burylo et al. (2011) also reported a high variability of RAR that was 
resulted from the heterogeneity in soil properties, the existence of natural obstacles, and also the 
plant interactions such as competition. Stokes et al. (2008) suggested that this variability was 
related to the interactions of genetic and environmental factors. Also, coarse roots can have a 
strong effect on RAR variability due to their larger cross sectional area. In general, most 
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researchers have reported a reduction in RAR values with increasing soil depths (Abernethy and 
Rutherfurd, 2000; Burylo et al., 2011; Abdi, 2014). This decrease has been attributed to the 
decrease of aeration, loss of nutrients storage, and the presence of more compacted layers 
(Bischetti et al., 2005; Abdi, 2014). Generally, the maximum amount of root density can be seen in 
the upper layers of soil near the surface (Burylo et al., 2011), therefore, fixation or reinforcement 
effects at these layers are expected to be more effective. In contrast with the previous researches, 
in this study, the maximum amount of root density was observed at the middle layers of soil 
instead of shallower layers. Genet et al. (2006) stated that a decrease in root biomass in shallower 
soil layers may be attributed to the deterioration in nutrient status and biological conditions. In 
our research area, a relatively high rate of wind erosion may cause unfavorable biological 
conditions in shallower soil layers and has caused a decrease in RAR in shallower layers. RAR 
values in uphill, downhill, leeward, and windward sides had no significant difference, except in 
two layers. The result may be due to the relatively smaller size of H. persicum height, crown and 
its roots compared to tree species that have been investigated in previous literature. The 
morphological characteristic of H. persicum can be considered as one of the advantages of this 
species, because as environmental stresses do not have a significant effect on modifying the root 
system, plantation spacing in different directions can be uniform and calculating the efficiency of 
afforestation for erosion control may be easier. 

It has been reported that the function of roots in the uphill, downhill, leeward, and windward 
sides is different because of various mechanical stresses (Sun et al., 2008). Therefore, in the 
current research, mechanical tests on tensile strength of roots were carried out in uphill, downhill, 
leeward, and windward samples. Based on previous studies, uphill and windward roots are 
reported to be stronger due to higher stress levels (Stokes, 2002; Abdi et al., 2010). However, no 
significant differences were observed in tensile strength between different treatments, just like 
RAR values. This may be related to the relatively small size and small crown of H. persicum that 
resulted in light weight and lower level of gravity and wind induced stresses. 

The relation between tensile strength and root diameter followed a negative power function, as 
reported in many previous studies (Abdi et al., 2010; Burylo et al., 2011; Vergani et al., 2012). A 
thin bark and a higher percentage of cellulose content in thinner roots can explain this 
phenomenon (Genet et al., 2005; De Baets et al., 2008). It should be noted that due to the different 
structure of the cellulose and lignin, cellulose is more resistant against tensile force (Genet et al., 
2005). 

Burylo et al. (2011) reported that the suggested range of a and £ coefficients of power function 
parameters for shrubby species was between 91.20>a>4.40 and —52.00>f>-1.75. According to 
the coefficients obtained from this study, all values are in the suggested range, except a 
coefficient for southward treatment that is lower than the suggested range (Table 1). H. persicum 
is an arid ecosystem species; however, the suggested ranges are based on European and 
Mediterranean ecosystems that may obtain the contradictory results. Therefore, the researchers 
need to modify the factors based on the species from other ecosystems. For any species, large 
amounts of a and small amounts of f, with respect to the relationship between strength and 
diameter, produce high tensile strength to prevent soil erosion, making it is a useful criteria for 
selecting species (Bischetti et al., 2005; Lateh et al., 2014) for restoration and erosion control. 
Although, information on arid species is scarce and future researches should aim to cover this 
gap. 

In this study, the model of Wu et al. (1979) was used to estimate increased soil cohesion due to 
the presence of H. persicum roots. The results showed that increased cohesion effect of H. 
persicum ranged from 0.001 to 1.400 kPa. The ranges for uphill and downhill sides were 
0.020-0.470, and 0.005-1.400 kPa, respectively. Furthermore, this ranges varied from 
0.001—0.590 KPa for northward and 0.001—0.650 kPa for southward. It is worth mentioning that 
most of the previous studies concerning tree species and information about shrubs, especially in 
arid ecosystems are scarce. The increased cohesion due to the presence of roots estimated in the 
current study is relatively lower than that in previous studies (e.g., Abernethy and Rutherfurd, 
2000; Mattia et al., 2005; Genet et al., 2008), but is significant higher than that of the sandy soil 
in the study area (i.e., 0.080 kPa). As sandy soil is susceptible to wind erosion due to a lack of 
cohesion among its particles, increasing its cohesion can have a great effect on its fixation rate 
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and the resistance against erosion. It is also worth mentioning that the strength of coarse-grain 
soil mainly depends on their internal friction angle or ø. As soil shear strength controls the rate of 
splash, rill and sheet erosion and also mass movement, increased cohesion due to the presence of 
roots in soil can prevent or decrease these catastrophic events. Previous researches have reported 
variations of the amount of increased cohesion under the influence of mechanical stresses (Stokes, 
2002; Abdi et al., 2010). Although no significant effect of mechanical stresses (i.e., wind and 
slope) was found on increased cohesion in the current study, uphill and downhill profiles had 
lower levels of increased cohesion than downhill profiles, while they have expected to bear higher 
level of stress from slope and wind, thus having a higher increased cohesion. This finding is in 
contrast with previous literature (e.g., Stokes, 2002) and showed that environmental stresses do 
not have significant effect on the level of increased cohesion of H. persicum. The small size and 
small crown of H. persicum that resulted in light weight and lower level of gravity and wind 
induced stresses, can explain this phenomenon. 


5 Conclusions 


In this study, the role of H. persicum in increasing soil cohesion and therefore increasing 
resistance against erosion was examined through the assessment of RAR and root strength (tensile 
strength). The RAR was low at shallower layers and the highest values were observed at the 
middle layers of soil. The results confirmed a positive relation between root diameter and tensile 
force. Also a negative relation was seen between root diameter and tensile strength. In addition, 
soil cohesion in different soil depths followed a similar trend with RAR. Moreover, slope and 
wind induced stresses had no significant effect on increased cohesion presented from H. persicum 
roots. As the soil of the study area was coarse grain soil with low cohesion, added cohesion due to 
the presence of roots, can have a significant effect on the resistance of soil to erosion. The results 
presented here expand the knowledge on the biotechnical properties of H. persicum species as an 
important species in arid land afforestation. It is worth mentioning that the lack of information 
about the behavior and properties of root systems of arid and semi-arid species has led to 
limitations in soil bioengineering applications in such ecosystems. These findings may be used in 
assessing the efficiency of bioengineering techniques, afforestation and restoration measures for 
erosion control in these ecosystems. 
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